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Excitation Spectra in a Heavy-Light Meson-Meson System * 

H. R. Fiebig a (LHPCollaboration) 
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A system of two static quarks, at fixed distances r, and two light quarks is studied on an anisotropic lattice. 
Excitations by operators emphasizing quark or gluon degrees of freedom are examined. The maximum entropy 
method is applied in the spectral analysis. These simulations ultimately aim at learning about mechanisms of 
hadronic interaction. 


1. INTRODUCTION 

While QCD has emerged as the foundation of 
nuclear physics the question of how the strong 
hadronic interaction arises from first principles 
still awaits explanation. Lattice hadron physics 
holds the most promise for answers. Mechanisms 
of hadronic interaction may already be studied in 
systems of hadrons with one heavy, static, quark 
each. In this case the relative distance between 
the hadrons is well defined, and studies of the 
interaction mechanism in terms of intuitive de¬ 
scriptions, like potentials, becomes possible. 

In principle, knowledge of the excited states of 
a two-hadron system allows the construction of 
an effective interaction. In practice, the extrac¬ 
tion of excited states from a lattice simulation is 
a very difficult problem. The canonical approach 
is to use plateau fits to effective mass functions. 
Bayesian inference 0, a pillar in many fields of 
science, has surprisingly been ignored by the lat¬ 
tice community until only very recently |2|. While 
the aim of this work is to learn about hadronic 
interaction mechanisms extraction of the spectral 
mass density via Bayesian methods emerges as an 
interesting subject in its own. 

2. HEAVY-LIGHT MESON-MESON OP¬ 
ERATORS 

Our intention is to test the strengths of differ¬ 
ent excitation mechanisms, in a two-meson sys¬ 
tem, as a function of the relative distance r. The 
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following operators are used 

^l(^) = fir,x-y (1) 

Q A (xt)j 5 q A (xt) Q) 

\ _ ) \ _ ) 

= Yl ^r,x-y (2) 
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Up. AA '{xt,yt) Wp. B , B (xt,yt) 

Q A {xt)^q B {xt) Q B ,(yt)~i 5 q A '(yt ). 
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The operator dq represents a product of two lo¬ 
cal pseudo scalar mesons. The fields Q,q de¬ 
note heavy and light quarks, respectively. Indices 
A, A', B , B' mean color, and Up are link variable 
products along straight paths P. In <h 2 color sing¬ 
lets, i.e. mesons, spread across a distance r. Exci¬ 
tations due to these nonlocal operators involve 
gluon exchange degrees of freedom. In contrast, 
excitations involve states where the mesons 
mostly communicate via qq exchange, at large r. 
Diagrammatic representations of diagonal corre¬ 
lator matrix elements Ca(t,to ) = ( < &!(£) < L(M) 
shown in Fig. [I] illustrate these points. We will 
loosely refer to <I>i and <I> 2 as emphasizing gluon 
and quark degrees of freedom (d.o.f.). 

3. TECHNICAL LATTICE ISSUES 

Light-quark propagators are computed from 
random Z 2 -source estimators which are non¬ 
zero on time slice to only. Heavy-quark propaga¬ 
tors are treated in the static approximation, i.e. 
as link products along straight temporal paths. 
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Figure 1. Diagrammatic impressions of the corre¬ 
lators Cn and C 22 representing quark and gluon 
exchange degrees of freedom, respectively. 


Gaussian smearing of the quark fields and APE 
fuzzing of the gauge field links is applied. 

Simulations are done on a L 3 x T = 10 3 x 30 
anisotropic lattice with a (bare) aspect ratio of 
a s /at = 3. We use the tadpole improved gauge 
field action of || with (3 = 2.4. This corresponds 
to a spatial lattice constant of a s ~ 0.25fm, aj 1 ~ 
800MeV. The (quenched) anisotropic Wilson 
fermion action is augmented with a clover ac¬ 
tion limited to spatial planes. Only spatial direc¬ 
tions are improved with renormalization factors 
u s = (D) 1/4 , while Ut = 1 in the time direction. 
The Wilson hopping parameter k = 0.0679 leads 
to the mass ratio m w /m p ~ 0.75. 

4. SPECTRAL DENSITY 

For an operator combination +t> 2$2 

the correlation function 

a (Mo) = ($J(*)M*o)) = (3) 

= £|<n|^(to)|0}| 2 e-^ <■*-**) (4) 

71^0 

has forward (uj n > 0) and backward (ui n < 0) go¬ 
ing contributions on a periodic lattice. We nor¬ 
malize those differently, defining 

exp T (w, t) = Q{u)e~ ut + e{-oj)e +u( - T ~ t) , (5) 


the spectral model 

/ +00 

dco p T (u)exp T (u>,t-t 0 ) (6) 

-OO 

is then expected to fit the lattice correla¬ 
tors. Specifically, the requirement F(pr\t,to) = 
C v (t, to) leads to a sum of discrete (5-peaks for the 
spectral density 

PT (w) = ^2 S(u> - LOn) |(n|<l„(fo)|0)] 2 x (7) 

riy£ 0 

[0(w n ) + 0(-w n )e _ “ nT ]. 

The peak-like signature should survive in the dis¬ 
cretized form of (0) 

k+ 

F(p\t,t 0 )~ ^ pk exp T (vk,t-to), (8) 

k—K- 

where A u> is a suitable interval, ojk = A uik, and 

p k = Aup T (uJk)- 

From the Bayesian point of view the spec¬ 
tral weights pk,k = K-...K+, as well as the 
values C v (t,to),t = 0...T — 1, of the correla¬ 
tion function are stochastic variables, p and C, 
subject to certain probability distribution func¬ 
tions. The conditional probability V[C <— p] 
of C, given p, is known as the likelihood func¬ 
tion. For a large number of measurements we 
have V[C <— p] oc e - * 2 / 2 , see 0 , where 

X 2 = i C v(ti,to) ~ F(pr\ti,t 0 )} x (9) 

ti ,t2 

r^ 1 (tl,t 2 ) [C v (t 2 ,to) - F(p T \t 2 ,t 0 )] 

is the usual \ 2 -distance between the data and the 
model, and r~ 1 (ti, f 2 ) denotes the inverse covari¬ 
ance matrix It accounts for the statistical 
dependence of correlator data between time sli¬ 
ces. Assuming minimal information (in the sense 
of 0 ) about p the Bayesian prior probability is 
P[p\ oc e aS where a is a parameter and 

S[p] = y l(pk- rn k - Pk In — (10) 
k V mk ' 

is the entropy relative to to, the latter being called 
the default model Jlj. By virtue of Bayes’ theorem 
|7] the so-called posterior probability is 

v[p <— C] = V[C v- p]P[p\/P[C] Ok e~ w W , (11) 
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for a given C, where W is defined as 

W[p]= X 2 /2-aS. (12) 

The most likely p is obtained by minimizing W[p]. 
The method of choice ||, so far, has been singu¬ 
lar value decomposition (SVD). Indeed, the func¬ 
tional W\p\ has a unique minimum jlj]. Stochas¬ 
tic methods, however, seem more in tune with the 
Bayesian perspective. In this spirit we have em¬ 
ployed simulated annealing to obtain the spectral 
density. The corresponding partition function is 

Z w = J[dp\e- 0wW[p] . (13) 

A Metropolis algorithm was used with local up¬ 
dates pk —> xpk where x is a random deviate with 
distribution pi{x) = xe~ x . For the annealing 
schedule a power law f3w(n ) = (Pi — Po) (n/N) 1 + 
Po turned out to be suitable. 

5. RESULTS 

We have explored the a-dependence of p using 
an artificial (mock) data set compiled from eigen¬ 
values of the 2x2 correlation matrix Cij(t,to). 
Cooling averages were taken during the last 200, 
of N = 2200, cooling steps during which the final 
Pw was kept constant. For tuning purposes it is 
useful to monitor the entropy load 

s/ w ■■= • (i4) 

{W) p w ^oo 

At constant default model m = 10~ 12 the entropy 
weight a was varied over 15 orders of magnitude. 
Two examples are shown in Fig. Remarkably 
the gross structure of the spectral density remains 
stable across 8 orders of magnitude of a. The 
second pair of frames in Fig. marks a kind of 
critical point as, beyond which the annealing ac¬ 
tion W becomes entropy dominated, the onset 
of ignoring the evidence V\C\. This results in a 
smoothing of the contours of the large-w peak, see 
Fig. I Even so, the fit to the correlation function 
is not greatly affected. We found that a good cri¬ 
terion for choosing the entropy weight parameter 
a is to tune S/W to about < 10 -2 . As is evident 
from Tab. |l|, near this value the linear relation¬ 
ship between S/W and a begins to break down. 



Figure 2. Entropy weight dependence. Correla¬ 
tors and spectral densities for an artificial data 
set, at relative distance r = 1 and constant de¬ 
fault model m = 10~ 12 , for a = 2 x 10 -5 and 
2 x 10+ 4 , respectively. 


To explore the default model dependence the 
correlation function Cn of the local meson-meson 
operator (!l|) was analyzed for a series of constant 
models m. The results for two values, m = 10” 12 
and m = 1.0, are shown in Fig. |^. Different en¬ 
tropy weight parameters, a = 100 and 1000, were 
chosen such that the entropy load S/W is about 
the same. Increasing m by twelve orders of mag¬ 
nitude changes p at the high-w end. In that re¬ 
gion spectral peaks are apparently not supported 
by the data. Thus the spectral density tends to¬ 
wards the default model. Again, absence of in¬ 
formation leads to a smoother curve. Note that 
the change in p for u> > 3 is numerically inconse¬ 
quential due to the presence of large peaks, also, 
the gross structure of the spectral function is sta- 
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Table 1 

The entropy load S/W for some entropy weight 
parameters a. Note the linear relationship in the 
region a < 2 x 10 +1 . 


a 

S/W 

2 x 10~ a 

1.4 x 10~ s 

2 x 1CT 2 

1.4 x 1CT 5 

2 x 10+ 1 

1.4 x 10~ 2 

2 x 10+ 4 

8.6 x lCT 1 


ble against changing to. On the other hand, the 
micro structure of like peaks in Fig. ^ differs no¬ 
ticeably. 

It is well known that the functional (|l2|) has a 
unique absolute minimum jlj. Any start config¬ 
uration used in the annealing process will result 
in a final configuration p that is merely close to 
the absolute minimum, since (3w < oo. Then, 
which features of p are independent of the start 
configuration? Figure |] shows two spectral func¬ 
tions obtained from a cold start at to = 10 -12 
and a random start. There is a main and a sec¬ 
ondary peak (1' and 1"). Repeated computation 
of p with varying random starts and annealing 
schedules shows that the two-peak feature is sta¬ 
ble, thus being a property of the data. We have 
looked at peak volumes and energies 

Z n = J dLup T (u) ) = |{n|4>„(f 0 )|0)| 2 (15) 

E n = Z n 1 J du>p T (u))uj (16) 

where 5 n = {w | u> € peak #n}. In Tab. || there 
are listed the peak volumes obtained from a sam¬ 
ple of 4 starts (1 cold, 3 random), and the aver¬ 
ages (ave) and standard deviations (sig). Assum¬ 
ing that the spectral density comprises two dis¬ 
tinct states, V and 1", the variances of Zy and 
Zy are « 10. Their combined volume Zyjy 
however has a variance less by one order of mag¬ 
nitude. Peak splits, and also spikes, can be caused 
by imperfect lattice data. In this light, interpret¬ 
ing Fig. 0 in terms of only one physical state, 
1 = l'U 1", as opposed to two, 1 = 1' and 2 = 1", 
is a possibility. More analysis work will be needed 
to decide this matter. 



Figure 3. Default model dependence. Correlators 
of the two-meson operator <3>i at relative distance 
r = 2. The spectral densities are obtained with 
constant default models of to = 10 -12 and to = 
1.0, respectively, at constant entropy load S/W. 


Amoung the features of p that are computable 
via spectral analysis are low-ca moment integrated 
quantities, like the peak volume Z n and energy 
E n , for some state n. These are interesting ob¬ 
servables because they may reveal aspects of the 
physics of hadronic interaction. For example, the 
relative strength of excitations due to operators 
$2 versus is measured by £ = ln(Z 2 /Zi), in a 
normalization independent way. We have plotted 
£ in Fig. [| as a function of the relative meson- 
meson distance r. As r decreases <b 2 -type excita¬ 
tions gain strength over dq-type excitations. In 
view of (0), d>, and Fig. [j] we may say that gluon 
exchange degrees of freedom move to dominate 
quark exchange d.o.f. at small r. We consider 
Fig. | a preliminary result, a more detailed anal- 
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Figure 4. Spectral densities from a cold and a 
random start, for at r = 2, also see Tab. [| 


Table 2 

A possible peak split. With 4 different annealing 
starts the peak volume © fluctuates between 
peaks 1' and 1". 


Start 

Z v 

Zvu i" 


1 

142.0 

170.9 

28.9 

2 

162.9 

167.9 

5.0 

3 

155.4 

170.2 

14.8 

4 

141.1 

171.0 

29.9 

ave 

150.3 

170.0 

19.7 

sig 

9.2 

1.3 

10.4 


ysis is in progress. Possibly, in the chiral limit, 
level crossing of £2 and e\ in Fig. || may occur [||. 

6. CONCLUSION 

Studies of hadronic interaction require ex¬ 
tracting excited states from lattice simulations. 
Bayesian inference proves a superior tool com¬ 
pared to standard plateau methods. Using an 
entropic Baysian prior integrated quantities, like 
the excitation strengths and energies of certain 
states, can be reliably computed. The excitations 
of a heavy-light meson-meson give insight into as¬ 
pects of strong interaction physics. At the current 
lattice parameters, quark degrees of freedom still 
dominate the physics of hadronic interaction. 



r [0.25fm| r [0.25fm] 

Figure 5. Normalization independent excitation 
strength ratio ( of gluon versus quark exchange 
degrees of freedom as a function of the relative 
meson-meson distance r (left). Also shown are 
the corresponding energy shifts £ 1,2 compared 
to the mass 2m of two noninteracting mesons 
(right). The error bars reflect annealing start 
standard deviations. 
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